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PACS 64 . 70 . Q — Theory and modeling of the glass transition 

PACS 45 . 70 . -n - Granular systems 

PACS 47.57.-s - Complex fluids and colloidal systems 

Abstract. - We add relaxation mechanisms that mimic the effect of temperature and non- 
equilibrium driving to the recently-proposed spiral model which jams at a critical density pc < 1. 
This enables us to explore unjamming by temperature or driving at pc < p < 1- We numerically 
calculate the relaxation time of the persistence function and its spatial heterogeneity. We disen- 
tangle the three different relaxation mechanisms responsible for unjamming when varying density, 
temperature, and driving strength, respectively. We show that the spatial scale of dynamic het- 
erogeneity depends on density much more strongly than on temperature and driving. 



Introduction. — Glass-forming-liquids, colloids, 
' emulsions, foams, and granular matter all develop sluggish 
and heterogeneous dynamics as they approach the onset 
of jamming. The slowing down of the dynamics in these 
■ systems with increasing density of the constituent parti- 
I cles, decreasing temperature, or decreasing the strength 
r of external driving forces is often summarized in the form 
of a jamming phase-diagram [1]. To date, most numerical 
studies of this diagram have focused on particulate mod- 
els such as sphere packings. From a theoretical point of 
r. view, however, simpler models are easier to understand. 
Here, we introduce a lattice model with a phase-diagram 
. (fig. [1^) that is similar to that of sphere packings, and use 
I it to study dynamic heterogeneities. 

Experiments on granular [2] and colloidal [3] systems 
show steady growth in dynamic heterogeneities as the re- 
laxation time increases with increasing density. In glass- 
forming liquids, however, the scale of heterogeneities re- 
mains modest even as the relaxation time increases by 
more than 10 orders of magnitude with decreasing temper- 
ature [4] . This difference may be due to the far greater dy- 
namic range measurable in glass- forming liquids ^3|. Our 
model, however, suggests that this difference signals a fun- 
damental distinction between jamming due to density as 
opposed to jamming by temperature or driving. 

At zero temperature and driving, sphere packings un- 
dergo a jamming transition with density [SHZ]- This tran- 
sition has a mixed nature in that the number of interacting 



neighbors per sphere jumps discontinuously from zero to 
the minimum number needed for mechanical stability, but 
there are diverging length scales [BHH]. Above the criti- 
cal density, the spheres may be unjammed either by rais- 
ing temperature above a glass transition into an equilib- 
rium state or by applying shear stress above a yield stress 
to drive the system into a homogeneous non-equilibrium 
steady-state. 

Model. — To construct a lattice model with a simi- 
lar phase-diagram, we start with a kinetically-constrained 
lattice-gas. In such models, occupied sites do not interact 
but the dynamical rules governing changes in occupation 
of a site depend on the occupation of neighboring sites. 
This leads to dynamics that are increasingly slow and het- 
erogeneous as the fraction of occupied sites increases, since 
increasingly larger regions are required to rearrange col- 
lectively [9l-[TT|. In such models, the fraction of occupied 
sites might either be interpreted as a density variable or 
(in spin versions of the model) controlled by temperature. 
Thus, temperature and density are equivalent variables 
in such models, and cannot be controlled independently 
of each other. In order to generalize to the case where 
temperature and density are independent control variables, 
we choose to associate the fraction of occupied sites with 
density, p, so that a conventional kinetically-constrained 
model corresponds to our zero-temperature model. Here 
we focus on the spiral model jT51fTl] , defined on the square 
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Fig. 1: (Color online) a) Jamming phase-diagram of density p, 
temperature T, and driving v. Jammed phase is the thick red 
line along the p-axis which terminates at pc. Blue arrows are 
the trajectories we investigate, b) The spiral model is defined 
by dividing the neighbors of each site on the square lattice 
into four pairs, labeled NE, SE, SW, and NW (see text), c) 
Divergence of relaxation time as p pc for T = and v = Q, 
with lattice size indicated in the legend. 



lattice such that the occupation of a site can only change 
if its (NE or SW) and (NW or SE) neighboring pairs (see 
fig-Hh) are completely empty. This model jams into a non- 
ergodic phase at pc ~ 0.705, which allows us to study its 
behavior at Pc < P < 1 with nonzero temperature or driv- 
ing. Such jamming-percolation models |15fll8) possess an 
additional important property: the jamming transition at 
Pc has a mixed nature; the fraction of stuck particles that 
cannot participate in rearrangements jumps discontinu- 
ously as in a first-order transition, while time and length 
scales diverge as in a second-order transition. For a re- 
cent investigation of the equilibrium and non-equilibrium 
dynamics of the spiral model, see [19] . 

In the original spiral model the number of particles 
was not conserved. Following |20j . we modify the zero- 
temperature stochastic dynamics so that instead of switch- 
ing sites between being occupied and vacant we move a 
particle to a neighboring site if the target site is vacant 
and if the kinetic constraint described in fig. [lb holds both 
before and after the move. We measure time in units of 
attempted moves per particle. 

We introduce temperature by softening the kinetic con- 
straints. Instead of preventing blocked moves, we allow 
them with probability exp(— l/T). Thus, for T = we 
recover the original model with rigid constraints. Note 
that the system still has no interactions and that energy 
is only associated with the virtual barrier the system has 
to cross in a kinetically-constrained move. For p > p^, the 
system is non-ergodic only at T = 0, since for arbitrarily 
low temperature, kinetically-constrained moves occur at a 
slow but nonzero rate, hence the system may eventually 
reach any configuration. 

We drive the system into a non-equilibrium steady-state 
by inducing a current of particles in one direction (for ex- 



ample, from left to right), as follows. In addition to the 
afore-mentioned moves in which particles can move into 
neighboring vacant sites subject to the soft (T > 0) or 
rigid (T = 0) kinetic constraint, we introduce a second 
type of move, in which a particle can move into the neigh- 
boring site to its right if it is vacant, irrespective of the 
kinetic constraints |21) . Such moves are attempted at rate 
/, and we characterize the driving strength by the aver- 
age flow velocity v = {1— p)f induced by them. Note that 
even for arbitrarily slow driving, the environment of any 
blocked particle will eventually change such that the par- 
ticle will no longer be blocked and can move even under 
the kinetically-constrained dynamics [24j . 

We could have chosen to soften the kinetic constraints in 
more complicated ways that couple density to temperature 
or driving to obtain more realistic results. The advantage 
of our implementation is that the interplay of density, tem- 
perature and driving appears in its purest, simplest form. 

We explored the model with rejection-free Monte-Carlo 
simulations, and established convergence of the results 
with system size by comparing systems of 80x80, 400x400, 
and 2000x2000 sites. To extract the relaxation time and 
the scale of dynamic heterogeneities, we calculate the per- 
sistence function, Pi{t), defined as the probability that 
particle i has not moved over a time interval t [25]. At 
long waiting times we find either exponential or stretched 
exponential decay, depending on density, of the particle- 
averaged persistence function, p{t) = i/N^^pi{t), where 
N is the number of particles in the system; we extract 
the relaxation time r from the condition p(r) = 1/e. 
We measure the dynamic heterogeneity in the standard 
way [2BH25] , by calculating the variance of the persistence 
function 

Xi{t) ^ N[{p\t)) - {pit))% (1) 

where ( ) denotes an average over different stochastic real- 
izations for the initial state and dynamics. This quantity 
has been shown to reflect the spatial extent of dynamic 
heterogeneities [51[30]. The idea behind this is that if re- 
laxation is heterogeneous, different regions in space, or 
alternatively different copies of the system, relax at differ- 
ent times, hence there is a large variance in p{t) between 
them. 

In the absence of temperature or driving, relaxation 
slows down with increasing density not only because each 
particle is more likely to be blocked by its neighbors but 
also because these neighbors are in turn blocked by their 
neighbors and so on. Hence moving a particle requires a 
collective motion of many other particles. In the spiral 
model, as p ^ Pc, the size of blocking clusters is expected 
to diverge as log(L) - {pc - p)"'' with p = 0.64 [l4l[T5] . 
Fig. [TJ; shows that we find log(T) ~ {pc — p)~'^- This is 
consistent with a scaling relation between length scales 
and time scales of the form r ~ L^, so that log(T) ^ 
log(L), irrespective of the scaling exponent z. 

Jamming Mechanisms. — The behavior of the re- 
laxation time as a function of temperature and density at 
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Fig. 2: (Color online) a) Relaxation time vs inverse tempera- 
ture for various densities, b) Relaxation rate due to the temper- 
ature mechanism at all densities may be collapsed to Arrhenius 
form by normalizing Rt by its T = cxd value given in the top 
inset, and scaling T by the effective barrier height A{p) given 
in the bottom inset. 



Fig. 3: (Color online) a) Relaxation time vs inverse driving 
for various densities (same legend as fig. b) Collapse of 
relaxation rate due to the driving mechanism for all densities 
by scaling _R„ by the « = oo value and scaling by («o/w)". 
R'^(p) and vo[p) are given in the inset and a[p) is given in 
fig. [St. 



zero driving is summarized in fig. [2^. As T 0, r diverges 
for p > while for p < pc '\i saturates to the finite T = 
value given in fig. [TJ;. The singularity at pc and T = af- 
fects the behavior of r(p, T) nearby. However, the overall 
relaxation rate we measure arises from a combination of 
two types of physical processes, which may be attributed 
to density and temperature separately. We demonstrate 
this by writing the relaxation rate as 



1 



Rp{p)+RT{p.T). 



(2) 



Here, Rp = 1/t{p,T — 0) is the relaxation rate due to 
the density mechanism, which represents processes sub- 
ject to the kinetic constraints in which the neighborhood 
of a particle changes so that a previously-blocked parti- 
cle can move. Thus, Rp = for p > pc- Similarly, Rt 
is the relaxation rate due to the temperature mechanism, 
representing the process in which a blocked particle moves 
either by directly overcoming the kinetic constraint or by 
becoming unblocked when one of its neighbors overcomes 



its kinetic constraint by thermal activation. Clearly, 
Rt 



t{p,T) Tip,T = 0) 



(3) 



Figure [2b shows that Rt has an Arrhenius dependence 
on temperature: Rt = R^ exp{—A{p)/T). Here, Rt{p) 
is the relaxation rate at T = oo, where the kinetic con- 
straint becomes irrelevant. For p > pc, the dominant pro- 
cess facilitated by thermal activation is motion of a single 
blocked particle moving at probability exp(— l/T), there- 
fore A 1. For p ^ Pc on the other hand, the density 
and temperature mechanisms are inherently coupled since 
moves are typically blocked by a large cluster of neighbors, 
and there are multiple moves that can lead to unblocking 
a single particle, hence A < 1. Such collective behavior 
only for p < pc and not for p > pc, which eventually leads 
to the non-trivial form of A{p), is another manifestation of 
the mixed (or one-sided) nature of the jamming transition. 

Figure [3^ shows the relaxation time vs driving for vari- 
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Fig. 4: (Color online) Ratio of the simulated relaxation time to 
the prediction of eq. ((5| vs the ratio of temperature to driving 
relaxation rates. 



ous densities at zero temperature. We now write 
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t{p,v) 



Rp{p) + Rv{p,v), 



(4) 



where Ry is the relaxation rate due to the driving mech- 
anism, in which the neighborhood of a blocked particle 
changes due to driving events that unblock it and enable 
it to move subject to the kinetic constraint. We obtain 
Ry by subtracting Rp = 1/t{p,v ~ 0) from 1/t{p,v) and 
show in fig. [Sb that Ry/ has the same dependence on 
{v[)/v)'^ for all p. Here, R^{p) describes the relaxation 
rate at infinite driving strength, vo{p) decreases mono- 
tonically with p, and a{p) behaves similarly to A{p) (see 
insets to figs-llb andlSla). Note that above pc and at densi- 
ties below Pc but away from the vicinity of the transition, 
a ~ 1. Thus, Ry varies linearly with the fiow velocity, v, 
at small v, as expected. At higher v, for all densities, Ry 
crosses over to a constant at v ~ vq because the neighbor- 
hood around a particle is completely randomized between 
attempts of diffusive moves. As a result, increasing the 
driving strength even more does not affect the relaxation 
rate at high flow rates. 

We now consider the interplay between temperature and 
driving. So far, we have identified the thermal relaxation 
rate, Rt, for v = 0, and the driving relaxation rate, Ry, 
for T = 0. For T > and v > 0, the simplest assumptions 
are that Rt does not depend on v, Ry does not depend 
on T, and the relaxation rates are additive, so that 



t{p,T,v) 



1 



Rpip) + RTip,T) + Ryip, v)' 



(5) 



Figure H] shows the ratio of the actual relaxation time mea- 
sured in simulations in which both T > and u > to this 
prediction. Obviously, when the values of Rt and Ry are 
very different, the smaller rate becomes irrelevant and the 
larger behaves as it behaves in the complete absence of the 
smaller. When Rt and Ry are comparable, deviations of 
around 20% are seen, indicating that the two mechanisms 
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Fig. 5: (Color online) Dynamic heterogeneity vs relaxation 
time along various paths to jamming, as indicated in the leg- 
end. 



are coupled and the relaxation is not given as a simple 
sum of independent relaxations. For p < pc, these devia- 
tions disappear for small values of Rt + Ry since then the 
density mechanism dominates. 

Dynamic Heterogeneity. — Differences between the 
relaxation mechanisms are clearly visible in spatial corre- 
lations of the dynamics [25^29) . For a given p, T, and v, 
X4{t) is maximal roughly when t = t with a value related 
to the typical number of particles that rearrange collec- 
tively [2j|30]. Figure [5] shows the maximal value of Xi 
vs T along different paths in the jamming phase-diagram. 
Solid circles denote the path of increasing density at T = 
and w = 0. As expected, along this path, the typical 
size of clusters rearranging collectively diverges as the sys- 
tem jams [15j . The data behaves as max(x4) oc with 
A w 0.5, in reasonable agreement with granular [2 and 
colloidal [3] experiments. 

When T decreases at fixed density, X4 does not diverge 
unless p — Pc- For p < pc both the density and tempera- 
ture relaxation mechanisms are at play. A blocked particle 
can move by waiting until its neighbors move to unblock 
it via the density mechanism. Additionally, temperature 
assists relaxation not only by allowing a blocked particle 
to overcome the kinetic constraint on its own, but also by 
allowing the neighbors of this blocked particle to overcome 
their kinetic constraints, thus releasing it and enabling it 
to move by a unblocked move. Overall, the temperature 
mechanism is less collective than the density mechanism. 
This is shown by the result in Fig. [S] that X4 is smaller 
(left-pointing triangles) at T > than along the T = 
path. As T decreases, the contribution of the temperature 
mechanism vanishes and the dynamics become dominated 
by the density mechanism and therefore become more col- 
lective with X4 increasing until it meets the T — 0, v — 
curve. For p > pc, relaxation occurs only via the tem- 
perature mechanism. Such relaxation involves primarily a 
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single blocked particle that waits a long time until it man- 
ages to move by a thermal move, and does not rely on the 
correlated dynamics of many particles. Here, X4 does not 
grow at all with decreasing temperature at p > (right- 
pointing triangles) , indicating that the typical spatial size 
of each rearrangement does not increase as the dynamics 
slow down. 

Experimental data for many glass-forming liquids shows 
that the scale of heterogeneities is also essentially constant 
with decreasing temperature at sufficiently low temper- 
atures. There, an initial rise with r in the number of 
molecules whose dynamics on the scale of r is correlated 
to a local enthalpy fluctuation is followed by a very slight 
increase or even saturation at large t [4j. 

When driving is lowered at fixed p < pc and T = (open 
squares in fig.[5|), both the density and driving mechanisms 
cause relaxation, and as for the thermal case, X4 grows un- 
til it meets the T = 0, w = curve. Interestingly, since the 
relative contribution of the density mechanism to the over- 
all relaxation is what determines the heterogeneity along 
the T -4- and w — ?> paths, it does not matter whether 
the non-collective rearrangement comes from temperature 
or from driving, and the curves for these two paths for 
p < Pc superimpose in fig. [S] For a path in which driving 
is lowered at p > pc and T — Q, the density mechanism 
is frozen out and relaxation is due to driving alone. As 
discussed earlier, the primary mechanism for relaxation 
in that case is that a blocked particle eventually becomes 
unblocked when its environment is changed by fiow. Like 
the temperature relaxation mechanism, this is a local pro- 
cess that does not become collective as the system jams. 
Fig.[5]shows that Xi is slightly larger along the u — ?> path 
than for the thermal case since more particles are involved 
in each rearrangement event, but the size of such events 
does not grow with increasing r (open diamonds), as for 
the temperature mechanism. 

Finally, we consider a trajectory aX, p > pc along which 
both temperature and driving are positive (stars in fig. [Sj . 
To maximize the interplay between the temperature and 
driving mechanisms we select T and v such that Rt — Rv ■ 
Since in this case the two non-collective mechanisms re- 
lated to temperature and driving govern the relaxation dy- 
namics, the spatial extent of dynamic heterogeneity sat- 
urates to a value which lies between that of the purely 
thermal trajectory (solid triangles) and that of the purely 
driven trajectory (open diamonds). 

Conclusions. — We studied a kinetically-constrained 
lattice-gas with a nontrivial jamming phase-diagram. Our 
model is substantially simpler than currently-used partic- 
ulate models with more realistic interactions, and is easier 
to study numerically so that a wide range of time scales 
(over ten decades) may easily be probed even in relatively 
large systems. 

Our model introduces three mechanisms by which den- 
sity fluctuations can relax, which we term the density, tem- 
perature and driving mechanisms. In particulate systems 



the same physical mechanisms come into play, but they 
are intermingled in a more complicated way. For exam- 
ple, real liquids, which are at high density, behave as if 
they have a lower effective density at high temperature. 
This is because increasing temperature in a particulate 
system increases the ability to open up free volume, and 
hence effectively decreases the density [3T]. It also de- 
creases the effective particle diameter, which effectively de- 
creases density [S^. Thus in real liquids, one would expect 
the effective density to increase as temperature decreases, 
thus slowing down relaxation due to the density mecha- 
nism. Once the density mechanism becomes too slow and 
is frozen out, the relaxation time should be dominated by 
the temperature mechanism, where small numbers of par- 
ticles overcome energy barriers via relatively uncorrelated 
rearrangements. Thus, it seems reasonable to expect that 
the sharp crossover from the density mechanism to the 
temperature one at pc in our model should be replaced 
in more realistic systems by a gradual crossover from a 
density-dominated regime at low densities or high tem- 
peratures, to a temperature-dominated one at high den- 
sity and low enough temperature. More generally, at suf- 
ficiently high densities and low temperatures or driving, 
the density mechanism should eventually freeze out, leav- 
ing only the temperature and driving mechanisms to relax 
the system. In that regime, our results suggest that the 
spatial scale of dynamic heterogeneities should be limited, 
as observed in glass- forming liquids {4]. Recent analyses 
of other lattice-based models find similar results [33ll34] . 

One corollary of this result is that one might expect the 
crossover from the density-mechanism-dominated regime 
at high temperature to the temperature-mechanism- 
dominated regime at low temperature to affect the form of 
the dynamics. It is known for hard-sphere systems, where 
only the density mechanism is at play, that the increase 
of relaxation time with decreasing temperature at fixed 
pressure is super- Arrhenius (fragile behavior) [3T]. On 
the other hand, the temperature mechanism should give 
rise to Arrhenius behavior (strong behavior). Thus, the 
crossover from a regime controlled by the density mech- 
anism to one controlled by the temperature mechanism 
might be accompanied by a fragile-to-strong crossover in 
the dynamics, although it is possible that one crossover is 
more abrupt than the other. 
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